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C.  1 


Abstract 

The  problem  of  the  plane  wave  oscillation  of  a  one- 
dlmenslonal  plasma  Is  considered.   Various  asymptotic 
expansions  of  the  solution  for  a  stable  plasma  are 
presented.   First,  an   expansion  is  given  for  short  time 
valid  for  arbitrary  wavelengths.   Second,  expansions  are 
given  for  long  wavelengths  valid  over  a  long  time  Interval. 
The  time  Interval  Is  sufficiently  long  to  observe  exponential 
damping  of  the  plasma  oscillations.   Whenever  the  oscillations 
are  strongly  damped  there  remains  an  Incoherent  part  of  the 
electric  field  that  slowly  and  Irregularly  tends  to  zero. 
Exponential  damping  Is  observable,  however,  only  in  the 
long  wavelength  limit. 

The  mixed  initial  value -boundary  value  problem  is  solved 
for  the  case  of  harmonic  time-dependent  data  given  on  a  wall. 
It  is  shown  that  for  stable  plasmas  the  electric  field 
approaches  a  limiting  form  in  time.   The  specular  reflection 
boundary  condition  is  solved  as  a  special  case.   The 
limiting  electric  field  is  characterized  by  an  analyticity 
condition  analogous  to  those  obtained  from  causality 
conditions  in  scattering  theory.   Finally,  asymptotic 
expansions  are  given  with  results  analogous  to  those 
described  above  where  the  small  parameter  analogous  to  wave 
number  is  the  deviation  of  the  frequency  of  oscillation 
from  the  plasma  frequency. 
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I.   Introduction 

One  of  the  simplest  problems  in  plasma  physics  is  that 
of  the  linearized  theory  of  the  plane  wave  oscillations  of 
a  plasma  composed  of  electrons  and  a  background  of  fixed 
ions  .   The  system  is  described  by  an  unperturbed  spatially 
homogeneous  distribution  function  of  electrons  f  (u),  the 
perturbed  electron  distribution  function  f(x,u,t)  and  the 
electric  field  E(x,t).   The  distribution  function  and  the 
electric  field  are  determined  by  the  linearized  Vlasov 
equation 

and  Poisson's  equation 


(2)  II  =  47rq  /f(x,u,t)du. 


The  initial  value  problem  for  a  plane  wave  solution  was  first 
solved  by  L.  Landau  in  19^5  ,  and  there  has  been  extensive 


1.  The  literature  on  this  problem  is  enormous  and  we  mention 
only  three  selected  references.   The  earliest  mention  of 
plasma  oscillations  is  in  a  paper  of  L.  Tonks  and  1. 
Langmuir,  Phys.  Rev.  33,  195  (l929).   L.  Landau,  J.  Phys. 
(U.S.S.R.  )  10,  25  (19^)  gave  the  first  satisfactory 
account  starting  from  the  linearized  Vlasov  equation. 

G.  Backus,  J.  Math.  Phys.  1,    I78  (1960)  has  examined 
this  problem  with  some  care  for  mathematical  problems. 

2.  See  1. 
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discussion  of  the  Interpretation  of  the  solution  ever  since. 
In  the  same  paper  the  boundary  value  protilem  for  the  system 
with  harmonic  time  dependence  of  r(x,u,t)  and  E(x,t)  was  also 
solved  for  a  special  f  (u)  and  with  specular  reflection  at  a 
wall.   In  this  paper  we  aim  to  clarify  some  questions  concerning 
both  problems  raised  by  us   and  by  others. 

In  his  discussion  of  the  Initial  value  problem  where 
f(x,u,t)  and  E(x,t)  have  spatial  dependence  e    ,  Landau 
considered  the  case  In  which  f  (u)  Is  a  Maxwell  distribution. 

0^  ' 

Under  the  assumption  that  f(x,u,o)  is  analytic  in  an 

appropriate  strip  in  velocity  space,  he  obtained  an  asymptotic 

expansion  of  the  electric  field  for  large  time.   This  expansion 

indicated  that  the  electric  field  oscillated  and  damped  in  time. 

For  small  k  (long  wavelength)  the  damping  was  very  small  and 

the  oscillation  frequency  was  near  the  characteristic  plasma 

frequency.   For  large  k  (short  wavelength)  the  damping  was  very 

strong  and  the  oscillation  frequency  was  far  from  the  plasma 

frequency.   There  were  many  later  "physical"  derivations  of 

the  oscillation  and  especially  the  damping  -Landau  damplng- 

which  obtained  analogous  results  for  cases  other  than  Maxwell 

distributions  and  analytic  initial  data.   However,  to  our 

knowledge  none  of  these  derivations  is  valid  nor  does  any  give 

5.   We  refer  to  three  of  our  earlier  papers  on  this  problem: 
I.   H.  Weitzner,  Phys.  Fluids  6,  1123  (1965),  II.   H. 
Weltzner,  Phys.  Fluids  -]_,    476  Tl964),  III.   H.  Weitzner, 
Comm.  Pure  Appl.  Math.  I8,  307  (1965).   We  refer  to 
these  papers  as  I,  II,  and  III. 
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a  full  description  of  the  behavior  of  the  electric  field.   We 
emphasized  In  an  earlier  paper  (l)  that  it  is  not  physically 
reasonable  to  consider  only  analytic  functions  for  initial 
data  and  that  without  this  assumption  the  electric  field  does 
not  have  an  asymptotic  expansion  of  the  form  described  by 
Landau.   We  showed  later  (ll),  however,  that  plasma  oscillations 
were  present  In  the  long  wavelength  limit.   We  were  unable  to 
find  a  complete  description  of  the  dispersion  of  the  plasma 
oscillations  in  time  and  we  said  that  Landau  damping  was  not 
observable. 

We  show  below  that  in  the  long  wavelength  limit,  but  only 
In  this  case,  it  is  possible  to  obtain  plasma  oscillations  and 
observable  Landau  damping  in  a  limited  time  interval. 
Asymptotically  in  time  no  oscillations  are  generally  present 
and  there  is  a  slow  irregular  dispersion  of  the  electric  field 
to  zero.   If  the  wavelength  is  not  long  then  our  previous 
results  apply,  and  after  an  initial  coherent  form  the  electric 
field  disperses  to  zero.   The  explicit  nature  of  the  dispersion 
of  the  electric  field  depends  extremely  sensitively  on  the 
detailed  properties  of  the  initial  data,  and  such  information 
is  not  physically  meaningful.   The  general  behavior  of  the 
solution  is  of  interest  in  fields  other  than  plasma  physics 
since  it  indicates  how  a  non-dissipative,  "collisionless" 
system  may,  in  appropriate  circumstances,  approximate  the 
damping  normally  associated  with  "collisional"  effects.   One 
certainly  expects  comparable  phenomena  in  other  problems 
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of  kinetic  theory  and  statistical  mechanics. 

For  the  problem  in  question  we  may  rewrite  the  system 
(1),  (2)  as 

where  (_J  is  a  linear  operator.   For  stable  plasmas  (we 
describe  the  concept  more  precisely  below)  (_J  has  only  a 
continuous  spectrum,  the  whole  Imaginary  axis.   But  in  the 
long  wavelength  case  there  is  an  approximate  point  eigenvalue 
on  or  near  the  imaginary  axis  which  induces  plasma  oscillations, 
Thus, in  this  limit,  the  problem  is  one  of  the  separation  of 
the  contributions  to  the  electric  field  from  the  neighborhood 
of  the  approximate  point  eigenvalue  and  the  remainder  of  the 
spectrum  and  then  the  estimation  of  their  behavior  in  time. 
We  are  now  able  to  extend  our  earlier  results  by  considering 
double  expansions  of  the  electric  field  in  terms  of  k  and  t 
In  order  to  find  a  more  complete  history  of  the  plasma 
oscillations.   For  a  short  interval  we  give  an  expansion  of 
the  electric  field   that  treats  all  parts  of  the  continuous 
spectrum  on  an  equal  basis.   If  we  consider  small  values  of 
k  in  this  expansion,  then  we  may  identify  terms  corresponding 
to  plasma  oscillations  and  contributions  from  the  remainder 
of  the  continuous  spectrum.   For  the  remainder  of  the 
expansions  we  take  k  small  and  t  as  specified.   For  t  of  the 
order  of  one,  or  not  too  large,  we  still  may  follow  the 
plasma  oscillations,  but  we  may  only  estimate  the  contribution 
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of  the  remainder  of  the  continuous  spectrum.   In  some  cases 
(but  not  for  Maxwell  distributions  or  other  with  Infinitely 
many  moments),  approximate  exponential  damping  In  time  may 
be  observable  In  this  time  Interval.   As  time  Increases,  but 
again  Is  not  too  large,  we  find  sharper  estimates  on  the 
contribution  of  most  of  the  continuous  spectrum  and  we  are 
able  to  observe  exponential  damping  of  the  plasma  oscillations, 
Finally,  for  sufficiently  large  time  the  coherent  plasma 
oscillation  has  dispersed  and  there  remains  only  a  slow 
Irregular  approach  to  equilibrium.   The  details  of  these 
expansions  are  given  In  the  second  section. 

The  boundary  value  problem  studied  by  Landau  Is  rather 
confused  In  that  even  for  the  special  case  considered  there 
are  no  Indications  of  proofs  or  ranges  of  validity  for  the 
expansions  that  exhibit  plasma  oscillations.   Rather  than 
consider  only  the  boundary  value  problem,  we  treat  the  mixed 
Initial  value-boundary  value  problem  and  then  approach  the 
limit  of  large  time.   We  take  the  prescribed  boundary  values 
to  have  the  time  dependence  e""""*^  .   In  the  special  case  of 
specular  reflection  we  obtain  the  well-known  result.   We 
choose  this  circuitous  procedure  to  select  the  one  causal 
solution  of  the  Infinitely  many  solutions  of  the  system  (l), 
(2)  with  time  dependence  e""^*^  .   We  show,  as  should  be 
expected,  that  the  solution  Is  uniquely  characterized  as 
the  one  which  possesses  an  analytic  continuation  In  the  full 
upper  half-plane  of  the  complex  variable  cd.   The  approach 
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to  equilibrium  applies  only  to  the  case  of  stable  plasmas. 
Finally  we  consider  asymptotic  expansions  of  the  solution 
corresponding  to  those  given  for  the  initial  value  problem. 
We  obtain  comparable  results  except  that  the  expansion  for 
small  X  no  longer  exhibits  plasma  oscillations  as  did  the 
one  previously  obtained  for  small  t.   We  give  the  solution 
to  the  initial  value-boundary  value  problem  and  discuss  the 
large  time  limit  in  Section  III.   In  Section  IV  we  examine 
the  expansions  that  show  the  plasma  oscillations. 


II.   The  Initial  Value  Problem 

If  we  assume  that  f(x,u,t)  and  E(x,t)  have  spatial 
dependence  e    ,  so  that  f  (x,u  ,  o  )  =  g(u  )e"^   ,  then  we  may 
solve  the  system  (l),  (2)  easily  by  taking  a  Laplace  transform 
in  time  and  eliminating  f(x,u,p)  in  favor  of  E(x,p).   After 
some  changes  of  variables  we  find 

00+1  a 

dzG(z)e-^^^^^ 


(3)       E(x,t)=  gqkA^e^^ 


(kA)2-p(z 


-00+ia 


2  - 
where  cd  =  (  '^^    )      is  the  plasma  frequency,  v  Is  the  thermal 

speed  used  to  non-dimenslonallze  the  distribution  functions  so 
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that 


f  (u)  -  ^  f  (^) 

O    '     V   O  V 


;(u)  =  ^  g(f ), 


n  is  the  electron  density  (chosen  so  that  /  f  (y)dy  =  l), 
A  =  —  Is  the  Debye  distance  and  "'" 


P 


(^)  G(z) 


r 

y-z    ' 


Imz  >  0 


\j 


—  oo 


5)         P(z)  = 


00 

^  Sf  (y)   H 

^  dy    ^  y-z  '      =^ 


The  contour  of  Integration  In  (3)  is  chosen  so  that  the 

integrand  Is  analytic  In  Imz  >  a.   We  show  below  that 

P(z)  =  Q(l/z  ),  and  hence  such  a  a  always  exists.   We  note 

that  the  original  equations  (l),  (2)  are  Invariant  under  a 

Gallllean  transformation  x  =  x-at^  t  =  t,  u  =  u-a.   Thus 

without  loss  of  generality  we  may  choose  the  reference  frame 

00 
In  which  /  yfQ(y)dy  =  0. 
—00 

If  P( z )  Is  real  and  positive  for  some  z  In  the  upper 
half -plane,  then  clearly  E(x,t)  grows  exponentially  In  time 
and  we  call  the  plasma  unstable.   For  our  purposes  It 
suffices  to  describe  a  plasma  as  stable  If  P(z)  Is  not  real 
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and  positive  In  Imz  >  0.   We  have  Included  the  real  axis  In 
the  definition  of  stability  and.  In  fact,  excluded  cases  in 
which  the  electric  field  Is  bounded  or  even  tends  to  zero, 
but  we  avoid  many  pathological  cases.   In  the  next  section 
we  admit  unstable  plasmas  in  considering  the  mixed  initial 
value-boundary  value  problem  before  discussing  the  approach 
to  equilibrium  for  stable  plasmas.   For  the  remainder  of 
this  section  we  consider  only  stable  cases.   The  Maxwell 
distribution  function  and  any  other  positive  distribution 
function  with  a  single  maximum  which  tends  to  zero  suf- 
ficiently rapidly  at  infinity  are  examples  of  stable  systems, 

h. 
although  many  others  exist. 

In  order  to  proceed  we  need  various  asymptotic 

expansions  of  P(z)  and  G(z)  for  large  |z|  in  Imz  >  0.   We 

omit  the  details  which  are  given  in  I  and  III  and  which  are 

quite  straightforward.   If  ^Ay)   has  N+4  moments  then  we 

readily  obtain 

N 


(6)        P(z)=^ 


"S:s-0(;fe) 


Where  3^  is  just  a  multiple  of  the  n   moment  of  f  (y)  and 


o 


the  expansion  is  valid  in  Imz  >  0.   If   f  (y)  has  infinitely 
many  moments  then  we  take  any  number  of  terms  in  the  series 


4.   0.  Penrose,  Phys.  Fluids  3,  258  (1960)  has  given  a  simple 
characterization  of  stable  plasmas. 


that  we  choose,  provided  |z|  Is  sufficiently  large.   For  real 
z  the  expansion  of  P(z)  has  no  imaginary  part.   The  presence 
of  an  Imaginary  part  is  essential  to  move  the  approximate 
zero  of  the  denominator  of  (3)  in  the  z  plane  from  the  real 
axis  to  the  lower  half-plane  and  thus  Induce  approximate 
exponential  damping  of  the  electric  field.   We  refer  to 
distribution  functions  with  infinitely  many  moments  as  type  I 
distribution  functions.   We  also  examine  a  second  type  of 
distribution  function  with  only  a  finite  number  of  moments  and 
with  the  asymptotic  expansion 

^-: rxy  -    Sgn(y)  ^/_     —^, 

in  which  case  F(z)  has  an  asymptotic  expansion  of  the  form 

Ivrp 
y  '^P^     A  +    1   \ 
"  ""   p    ^™P  U^    (-z)pJ  ' 


(if  p  is  an  Integer  then  the  term  In  the  second  summation  is 
replaced  by  -  ^  j(logz  -  lir)/ z^   +  ( logz  )/( -z  )P  J  ).   We  note 
the  presence  of  an  imaginary  part  in  these  series  for  real  z, 
Of  course  G(z)  has  similar  expansions  which  for  convenience 
we  take  of  the  form 


N 
a 


(8)  G(z)=zi^-Ori 


,   ..   ^,   N+l 
n=l  z      \  z 


Finally,  provided  ^^(y)  and  g(y)  have  many  well-behaved 
derivatives,  P( z )  and  G(z)  have  derivatives  on  the  real 
axis  and  we  may  use  Taylor  series  expansions  about  any 
point  on  the  real  axis  and  the  expansions  are  in  power  of 
(z-z  ).   Thus,  up  to  specified  errors  we  may  replace  P(z) 
and  G(z),  which  are  analytic  only  in  the  upper  half -plane, 
by  approximations  which  possess  continuations  into  the 
lower  half -plane,  and  then  attempt  to  evaluate  the  integrals. 

We  now  proceed  to  the  simplest  of  the  expansions.   We 
take  kvt  small  and  do  not  specify  kA.   We  choose  a  =  (kvt)   , 
so  that  z  is  uniformly  large  on  the  contour  of  integration  and 
we  may  expand  P(z)  and  G(z)  in  asymptotic  series.   We  write 
the  solution  for  the  explicit  case  of  three  terms  in  the 
expansion  of  P  and  G  and  the  general  error  estimate.   We 
find 

P   jw     \  (cXnZ   +a    )cos(kvtz    )-ia   z      sin(kvtz    ) 
(9)    E(x,t)  =  47rqkA2e^^    \       ^   ^     ^.  ^ ^^ ^- 

(a-,  z^-a^  )cosh(kvtz    )-ia_z    sinhfkvtz    ) 
1    c      3  ^ c  ^         2    c ^ c_ 

/-^  /  (kvt)^+^      f  ,kvt^2 

-    10    - 


where 


(10a)  zl  =   — ^  (l  +  /l+^p  (kA)' 


(10b)  z^  =  — ^-^  1^  l+4p2(kA)2  -  1  ^  , 


To  obtain  the  above  error  estimate  we  take  the  asymptotic 
expansions  of  P(z)  and  G(z)  up  to  and  Including  the  terms 
0(  l/z  ),  replace  the  denominator  In  the  error  term  by  Its 
minimum  and.  Integrate  the  rest  of  the  error  term.   Thus,  for 
any  kA,  If  kvt  Is  sufficiently  small  and  If  we  take  enough 
terms  In  the  series,  (9)  gives  an  explicit  expansion  of  the 
electric  field. 

Neither  the  original  expression  for  the  electric  field, 
(5),  nor  the  expansion  Is  regular  for  (kA)  small,  and  we 
examine  this  limit  with  much  greater  care  below.   For  (kA) 
of  order  one  or  larger  we  believe  that  the  only  general 
properties  available  are  (9)  and  the  result  obtained  by 
setting  a  =  0  In  (3)  and  repeated  Integration  by  parts: 


(11)   E(x,t)  =  WqkA^e^^''  Q 


(kvt)"  mln|(kA)2-P(z)i"+^ 
z  real 


For  kA  of  order  one  or  larger  we  have  the  expansion  for  (kvt) 
small  and  the  estimate  (ll)  for  kvt  large.   In  these  cases 
(9)  does  not  exhibit  oscillations  observable  over  several 
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periods  and  hence  does  not  represent  a  true  oscillation. 
The  estimate  (ll)  for  kvt  large  exhibits  the  dispersion  of 
the  coherent  form  of  the  electric  field  (9)- 

In  the  remainder  of  this  section  we  take  (kA)  small. 
In  this  case  there  are  two  roots  of  the  asymptotic  expansion 
of  (kA)^  -  F(z)  that  are  Q(l/kA)  and  all  the  other  roots 
are  0(l)  ^^  ^'^'      Clearly  tz   above  represent  the  roots 
(2)(l/kA)  and  ±z   are  the  roots  (Ji^)    in  kA.   We  note  that 
kvtz  -   Go(k)t,  where 

(12)  a:2(k)  =  CD^  +  ^^(kv)^ 

so  that  the  first  term  In  (9)  corresponds  to  plasma  oscillations, 

Further,  ±z   correspond  to  the  approximate  zeros  of  (kA)  -F(z), 

or  the  approximate  point  eigenvalues  of  the  system.   Also  ±z 

correspond   to  the  lowest  order  contributions  of  the  part  of 

the  continuous  spectrum  away  from  the  approximate  point 

eigenvalues.   We  denote  by  E  (x,t)  those  terms  In  (9)  that 

correspond  to  the  roots  +z  . 

P 

If  we  consider  type  I  distribution  functions  then  the 
expansion  of  (kA)   -  F(z)  to  any  order  Is  a  polynomial  In  l/z 
with  real  coefficients  and  hence  z   can  never  aqulre  an 
Imaginary  part.   If  we  consider  type  II  distribution  functions, 
then  the  expansion  to  a  sufficiently  high  order  will  have 


-kvt5 
v^,u;  ±t>  mu±b±pj.iea  oy  e 

where  to  lowest  order  In  kA 


an  Imaginary  part  so  that  E  (x,t)  Is  multiplied  by  e" 
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(13)  5  =  -  ^-^     ° 

(kA) 


2  /,.^n3  'Su 


u 


=  (kA)-l 


The  above  expression  Is  the  usual  form  of  the  damping  present, 
and  our  primary  task  in  this  section  is  to  show  that  for  (kX) 
small  there  is  a  time  Interval  In  which  such  damping  is 
generally  observable.   Surely  for  kvt'and  kA  both  small  no 
such  damping  is  significant. 

We  now  proceed  to  the  expansion  of  (j)  based  on  the 
assumption  that  kA  Is  small.   Initially  we  take  kvt  no  more 
than  of  order  one  or  somewhat  larger.   We  deform  the  contour 
of  integration  to  that  shown  in  Pig.  1,    where  we  also  define 
points   -Z  .   We  treat  separately  the  Integrals  from  -Z  to 
+Z   and  the  Integrals  from  ±Z   to  infinity.   In  principle 
we  should  either  use  molllfiers  or  carry  the  end  point  terms 
in  order  to  prove  that  the  end  point  terms  do  not  contribute 
to  the  results  below.   Either  technique  easily  shows  that  no 
end  point  terms  occur  and  we  omit  the  details.   Since  (kA) 
is  small  we  may  take  Z..  large  and  use  the  asymptotic  expansions 
for  |z|  >  Z  .   In  the  Integral  from  -Z^  to  +Z  we  Integrate 
by  parts  n  times,  replace  the  denominator  by  its  minimum  and 
estimate  the  resulting  Integral.   On  the  other  sections  of 
the  contour  we  use  the  expansions  of  F(z)  and  G(z),  picking 
up  an  error  term,  deform  the  contour  into  the  lower  half- 
plane,  picking  up  a  contribution  from  the  poles  at  iz  ,  and 
we  then  integrate  by  parts  n  times.   We  find  easily  that 
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(14)     E(x,t)  =  Ep(x,t)  +  ^—^-j  Q^__-_^^kqA 


(kvta)     /  ^ 
/v,  ,  v2  \2  0/"l^ 


log  Z 


^kA)2-i^y-^(kvtz^)" 


"i 


where  E  (x,t)  corresponds  to  the  oscillating  part  of  (9)  and 
has  the  damping  factor  exp(-5  kvt )  for  type  II  functions. 

We  may  use  (l4)  to  consider  many  special  cases  by 
appropriate  choices  of  a  and  Z  .   For  kvt  of  order   one  or 
smaller  we  take  o  =   1,    Z     ==  (kA)~  ,  1>  r>  0  and  we  readily 
verify  that  E  (x,t)  is  the  largest  term  in  (l4).   While  the 
second  error  term,  which  comes  from  the  continuous  spectrum 
not  near  z  ,  is  absolutely  large  it  is  small  relative  to 
E  (x,t).   At  a  later  time  this  term  is  absolutely  small  as 
well.   We  note  that  we  no  longer  have  an  explicit  approximation 
for  the  contribution  of  the  major  part  of  the  continuous 
spectrum,  and  we  only  have  the  estimate  given.   We  might  use 

(14)  to  follow  the  transition  from  kvt  of  order  one  to  kvt 
large,  but  we  omit  the  intermediate  values  and  take  kvt  large 
but  not  too  large.   We  choose  a  -    (kvt)    and  the  previous 
value  for  Z^ .   If  (kvt)  is  in  the  interval 

(15)  (kA)   ^  »  (kvt)  »  (kA)   ^, 

then  the  error  estimate  is  absolutely  small  and  the  plasma 
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oscillation  contribution  is  dominant,  provided  it  has  not 
damped  too  much  for  type  II  distribution  functions.   Again 
type  I  functions  exhibit  no  damping.   For  type  II  functions 
damping  is  observable  since  bf   /+v  ^vl/|v|^  and  kvtb  may  be 
large  and  yet  (15)  may  hold.   Even  in  this  time  interval  the 
amplitude  of  the  plasma  oscillations  may  become  so  small  that 
the  error  terms  dominate  since  E  (x,t)  damps  exponentially  fast 
and  the  error  terms  only  damp  algebraically. 

We  now  proceed  to  consider  the  limit  of  (kvt)  very  large 
and  find  damping  generally.   For  type  II  functions  this 
expansions  and  the  preceding  one  may  overlap.   We  move  the 
entire  contour  of  integration  onto  the  real  axis  and  we  define 
intervals  of  length  2t]  centered  on  C,^  where  Re  F(C,  )  =  (kA)  , 
C   >  0  >  C  }    and  C   are  large.   We  treat  these  intervals  and 
the  remainder  of  the  real  axis  separately.   Prom  the  parts  of 
the  Integral  away  from  C ,  we  find  after  n  integrations  by 
parts  (again  we  drop  end  point  terms  which  cancel  out)  that 
the  contribution  to  the  electric  field  is 


E,  =  qkA^  O  ^ 


(kvt)'^liE^  +  Ti(kA^)|"+l 


where  e  =  Im  F(C, ),  and  we  have  used  the  Taylor  series  for 
F(z)  near  C   to  estimate  the  smallest  value  of  the  denominator. 
We  may  approximate  the  remainder  of  the  integral  as 
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E   =  2qkA^   \/ 


+  kqA^OP^nr^(kA)^ 


L 


■le+CkA)-^!] 


3„i2 


We  now  deform  the  path  of  Integration  into  the  lower  half-plane 
and  pick  up  the  contributions  from  the  poles  near  ^  .   We  note 
that  there  is  no  difficulty  with  type  I  distribution  functions, 
but  for  type  II  distribution  functions  we  may  not  take  as  many 

terms  in  the  Taylor  series  as  we  please.   For  in  the  latter 

(n  ) 
case  e  /F^   (C, )  will  be  large  for  n  sufficiently  large, 

corresponding  to  a  large  shift  in  the  zero  of  the  denominator 

from  ^  .   Hence  in  this  case  we  must  restrict  n.   Finally  we 

integrate  by  parts  many  times  and  duplicate  the  previous 

error  estimate  to  find 


(16)    E(x,t)  =  E  (x,t)  +  r){ ± (kqA^ 

^  '  •    -.n.  2,  2/,  .^6^n+l  * 


(kvt)^He^+^^(kAr)- 


+  (kqA^)  Q(kA,)N+l(kA)g 
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where  we  now  Include  the  exponential  damping  with  decrement 
6  of  (13).   If  we  take   k:vt(kA)  r]  large  and  'ri(kA)   small  but 
much  larger  than  £  (as  we  may  easily  do),  then  the  error  terms 
are  small  and  the  net  damping  kvt5  =  (kvt  )£(kA)~'^  may  be   0(l) 
and  hence  observable.   Of  course  when  (kvt5)  becomes  very 
large  the  plasma  oscillations  have  damped  out  and  all  that 
Is  left  is  the  slow  dispersion  of  the  electric  field  to  zero 
for  which  a  best  estimate  is  Qjl/(kvte)  /   .   Thus  in  a 
limited  time  interval  and  when  kA  is  small,  exponential 
damping  of  plasma   oscillations  is  observable. 


III.   The  Mixed  Initial  Value-Boundary  Value  Problem 

In  this  section  we  consider  the  original  set  of  equations 
in  the  half -space  x  >  0.   At  the  wall  (x=0)  we  prescribe  the 
electric  field  and  the  distribution  functions  of  particles 
which  leave  the  wall.   To  complete  the  boundary  conditions 
we  require  that  the  disturbances  all  be  bounded  for  x  suf- 
ficiently large  for  any  fixed  t.   The  initial  conditions  is 
just  that  imposed  for  the  infinite  problem;  we  give  the 
distribution  function  of  particles  for  all  x.   These 
conditions  are  reasonable  physically,  and  we  see  below  that 
they  lead  to  a  well-posed  mathematical  problem  We  carry 
out  the  solution  explicitly  only  for  the  case  of  harmonic 
time  dependence  of  the  data  given  on  the  wall,  although  the 
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methods  obviously  apply  more  generally.   We  show  that  for  a 
stable  plasma  the  electric  field  has  a  limit  as  t  tends  to 
infinity  and  that  the  limiting  electric  field  is  uniquely 
characterized  by  the  prescribed  boundary  data  and  an 
analytlcity  condition  analogous  to  that  obtained  from 
causality  conditions  in  scattering  theory.   At  the  end  of 
this  section  we  obtain  the  solution  for  a  specular  reflection 
boundary  condition  as  a  special  case. 

We  take  the  initial  data  and  conditions  on  the  boundary  as 

f(x,v,0)  =   g(x,v) 

rh+(v)e-^^^       ,  V  >  0 
(17)     i      f(0,v,t)  -) 

h"(v,t)  ,v  <  0 


V_ 


E(0,t)  =   E   e"^"'^ 


where  all  the  functions  are  given  except  h~(v,t),  which  is 
determined  by  the  condition  that  E(x,t)  and  f(x,v,t)  be 
bounded  as  x  tends  to  infinity  at  any  fixed  time.   We  perform 
Laplace  transforms  of  the  system  (l),  (2),  both  in  time  and 
space  and  let  p  represent  the  transform  variable  corresponding 
to  t  and  s  represent  the  transform  corresponding  to  x.   We 
expect  the  double  transforms  E(s.,p)  and  f(sjV,p)  to  exist  and 
further  for  Re  p  sufficiently  large  we  expect  E(s,p)  and 
f(s,v,p)  to  be  analytic  in  Re  s  >  0  and  have  at  most  simple 
poles  on  Re  s  =  0.   We  do  not  consider  the  analytlcity 
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condition  on  f(s,v,p)  further  below,  for  an  elementary 
computation  shows  that  If  E(s,p)  Is  analytic  In  Re  s  >  0 
then  f(s,v,p)  Is  as  well.   The  transforms  of  the  system  (l), 
(2)  are 

r 


Sf  (o) 
(p+sv)f(s,v,p)  +  ^  E(s,p)    g^   =  g(s,v)  + 


vh  ( V ) 
p+lo) 


V  >  0 


vh  (v,p),   V  <  0 


00 

r 


v_ 


sE(s,p)  -  ^Trq 


f  ( s,v,p)dv  = 


E 


o 


p+ico 


<J 


On  elimination  of  f(s,v,p)  from  the  equations  and  with  the 
definitions 


f 


(18)        F(z)  = 


dv  \    o 


Sfjv) 


W+z '^    dv 


J 


(19)        G(z,s) 


^^^\    (^^)  g(s,v) 


(20) 


r 

H(z)  =  4Trq   (^)  vh+(v) 


v+z 


J 


o 


(21)        K(z,p)  =  47rq 


ip-)   vh-(v,p) 

v+z  \         i  ir-  J 


—  00 
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we  find 


(22)      E(s,p)=— ^±i^ 


sE  +H(p/s) 

°       -  +  K(p/s,  p)  +  G(p/s,  s) 


s^  +  cOp  P(p/s) 

Clearly  F(z)  has  an  expansion  of  the  form  (?). 

We  see  from  (l8)  -  (2l)  that  P(z)  and  G(z)  are  defined 

for  all  values  In  the  complex  z  plane  except  for  real  z  and 

that  P(z)  and  G(z)  have  cuts  on  the  real  axis.   Further  H(z) 

(k(z,p))  has  a  cut  on  the  negative  (positive)  real  axis,  and 

both  functions  are  bounded  at  the  origin  and  tend  to  zero  at 

infinity.   As  (22)  stands,  E(s,p)  is  discontinuous  across 

p/s  =  z  real  and  positive,  and  E(s,p)  may  be  singular  at  any 

of  the  zeros  of  the  denominator  of  (22).   We  must  choose 

h~(v,p)  so  that  in  Re  p  >  k.  Re  s  >  0,  E(s,p)  is  analytic  in 

both  variables.   We  may  write  the  denominator  of  (22)  as 

p 
•^  (p  +  CD   (p/s )  P(p/s ) ).   Since  |z  P(z)|  is  bounded  by  some 

constant  it  follows  that  for  Rep  sufficiently  large  E(s,p) 

has  no  poles.   We  ultimately  extend  our  expression  for  E(s,p) 

to  all  p  in  Re  p  >  0  for  stable  plasmas.   In  order  to 

suppress  the  cut  on  z  =  p/s  real  and  positive  it  suffices  to 

require  that  E(-^,p)  be  continuous  across  z  real  and  positive. 

We  may  restate  this  condition  as  a  Hil^ert  problem  for  K(z,p) 
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(23)     K+(z,p)  =  )     I     I  K-(z,p) 

p  +0)  z'^P~(z) 


+ 


E^  p/z  +  H(z)M  /p^^gzV(z) 


p'^+a3^z'=^P  (z)  )        ^ 

z  real  and  positive, 

where  the  ±  refers  to  the  limit  taken  from  above  or  below  the 
real  axls« 

To  solve  (23)  we  start  by  defining 


(24)         r(2)  = 


27rl 
•J 


00 

^  ^,        /p^+aj^t^P^(t) 

t-z^   "^^    2,  2,2^-/,  N 
P  ^t  P  (t) 


o 

and  In  (24)  we  select  the  branch  of  the  logarithm  that  tends, 
to  zero  as  t  tends  to  Infinity,   Again  for  Rep  sufficiently 
large  the  logarithm  tends  to  zero  as  t  tends  to  zero  as  well. 
Clearly  e'  ^'^'^^    satisfies  the  homogeneous  form  of  (23)  and 
the  standard  arguments  show  that  no  solution  of  the  homogeneous 
equation  Is  bounded  at  the  origin,  tends  to '-zero  at  Infinity 
and 'Is  analytic  except  for  a  cut  on  the  positive  real  axis. 
Thus  the  unique  acceptable  solution  of  (23)  is 
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(25) 


K(z,p)  = 


,r(^.p) 

27rl 


r  ,,,(Apa)  e +H(t) 


■J 


dt 
t-z 


p  +  Ico 


/^-r-(t,p) .  3-rnt,p)Y  G-(t,  pA)e-n"(t.p) 


-  G+(t,  p/tje-r'^'t.P) 


We  now  wish  to  modify  the  contours  of  Integration  In  (24)  and 
(25)  so  as  to  exhibit  explicitly  the  desired  analytlclty 
properties  of  E(s,p)  and  to  obtain  representations  appropriate 
to  study  the  approach  to  equilibrium. 

We  rewrite  (24)  as  two  Integrals  by  the  identity 

■p^4<i)^z^P^(z) 


log 


p^4<u^z^F^(z) 

p^-Ki^^^Sp-^^) 


=  log 


p 


2   2 
p  -Ro 


^p^+co^z^F  (z) 
-  log  I § — 5 ;  > 


2,  2 

p  +C0 

^   p 


and  each  Integral  converges  separately.   We  take  the  first 
(second)  Integral  and  deform  the  contour  of  integration  to  a 
line  perpendicular  to  p  in  the  upper  (lower)  half  plane  and 
introduce  the  variable   t  =  ipt  to  find 


(26a)    r(z,p)  =  2^ 


00 


dt 


—00 


(t    -z/lp) 


p2(l-co^(t    )%(lpt 
log  I  f-^ 

p  ^ 


+   log 


r        O  0        0 

p  +0)  z''P(z) 

2      2 
p    +C0 
P 


Ti(Re   p/z). 
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or 


(26b)      r(2,p)  -  r{z,p)  +  ic 


fp^+co^z%(z  ) 
2,  2 

P 


"Ti(Re  ^), 


where  r](x)  Is  the  Heavislde  step  functlono   We  perform  a 
similar  manipulation  on  all  the  terms  in  (25)  with  the 
exception  of  the  term  proportional  to  E  for  which  we  use 
the  contour  of  integration  of  Pig.  2  and  we  obtain 


(27)   K(z,p) 


Q^f ._.  .-r(z.p) 


+ 


27ri 


27ripE  e 


[-dvi.    e 


-r(o,p) 


pE 


■^  o    ,  ti[  z  J  ,  „r        P > 

Z  (  p+lCO  j     p+160     "•   '  Z  ' 


z  ( p+io) } 


u 


r  dt'e-r(iPt  .P)  ^H(l£)  ^  ,^,p,._  _iJ 


t  -  z/lp 


o 


In  (27)  we  have  taken  Re(p/z)  positive.   From  the  definition 
of  E(s,p),  the  relation  p  =  zs  we  infer 

(0,p) 
,  p+ico") 


^r(p/s,    p)-log(p2+cD|)    (27ri   p^E     e- 
(28)      E(s,p)   =^ > 


Siri 


"    J     V       (st'+l)  J\P*^'^  ^'^  It'Jj 

—00 

In  (28)  we  have  exhibited  E(s,p)  as  an  analytic  function  of 
s  and  p  in  Re  s  >  0,  Re  p  >  k,  so  that  (28)  represents  an 
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explicit  solution  to  the  problem.   Before  proceeding  further 
we  note  an  Identity  derivable  from  (27)  by  taking  the  leading 
term  in  the  asymptotic  expansion  for  large  z 


-1  °°      4-  °       - 

(29)  ^^(P^o^^^'^/  vh^(v)dv)   +   47rq   /   vh    (v,p)dv 

■^  O  -00 


PE     e-n(°.P) 

+   47rq/dxdvg(x,v)   -  — °   p^^^ 


00 


2Tr 


^dt'e-n(iPt',p)fHil£t:i^  G  Apt',      1 

\     p+ico  \    ^      'it 


We  might  use  (28)  and  (29)  to  study  lim  E(x,t)  Idj   obtaining 

X  -^  00 

the  residue  of  (28)  at  s  =  0,  but  the  results  are  none  too 
illuminating.   In  general  this  limit  is  not  zero  since  Eq. 
(2)  and  the  boundary  conditions  show  that  even  at  the 
initial  time  the  electric  field  is  already  not  zero  at 
infinity.   We  might  easily  verify,  however,  that  even  for:  stable 
plasmas   lim   lim  E(x,t)  ^      lim    lim  E(x,t). 

t— >  00  X— >oo  X — ^00   t— ^  00 

We  now  restrict  ourselves  to  stable  plasmas  and  we 
examine  the  approach  to  equilibrium.   We  see  that  if  P(z) 
is  never  real  and  positive,  the  condition  for  a  stable  plasma, 
then  (26)  provides  an  analytic  continuation  of  [''{p/s,    p) 
into  the  half-plane  Re  p  >  0,  p  =  ±  ico  excepted,  provided 
Re  s  >  0.   It  also  follows  that  (26)  provides  an  analytic 
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continuation  of  n(ipt  ,p),  t   real,  In  the  same  domain  in  the 

p  plane.   Thus  E(s,p)  Is  now  well  defined  In  Re  p  >  0, 

Re  s  >  0,  p  7^  -Ico,  +lcD  .   Finally  we  note  that  the  combinations 


(r(p/s,  p)  -  log(ajp+p^)  -  r(0,p))and  (TCp/s,  p)  -  log(p2+a)p 


P 
-  n  (ipt '  ,p))  are  well  defined  on  the  Imaginary  p  axis  at  all 

points,  and  these  combinations  are  dlff erentiable  with  respect 

to  p  except  at  p  =  tlco  .   Thus  E(s,p)  Is  defined  in  R  p  >  0 

except  at  p  =  -Ico,  where  It  has  a  pole,  and  the  standard 

arguments  show  that 


(30)      llm    (E(s,    t)ei^t)   =    -  ^  e  ^^  ^ "  ^    '    -la))-log(a32-c.2 ) 

t     -^    00 


-r(0,-la3)      _1 


00  , 

'"jj-'tt/     ^'>     -P(a3t     ,-103) 
dt   H(cDt    Je    '     ^         •*         ' 


■N 


2Trl 
J 


(st  +1) 


We  study  (30)  in  greater  detail  in  the  next  section,  we  note 
that  we  may  obtain  the  limiting  form  of  (30 )  for  large  x 
exactly  as  we  found  the  form  for  large  to   The  dominant 
singularity  is  at  s  =  0,  hence 


llm         llm   E(x,t)e^^^) 

X    -^  oo     t    -^   00 


-r(o,-icD)^  _^ 


00 

r 


+ 


ZIT 


dt'e-r'(^t',-ia3)j^(^^') 


A 


>  } 


u 


which  by   (29)   becomes 
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E      -   ^ 

O  0)1 

(31)        11m        11m   (E(x,t)e^'"^)   =  ^"^ 


vh   (v,t=oojdv 


2 

X    -^   00    t    ^    00  CD 

1     -    — ^ 
^  2 

CD 


Thus  with  no  current  flow  at  the  wall,  e.g.  for  a  specular 

reflection  boundary  condition,  the  limiting  electric  field 

is  E  (l  -  CD  /cD  )   ,  and  otherwise  there  is  an  additional 
o^     p'   ^   ■* 

contribution  from  the  current  flow  at  the  wall.   Eq.  (31 )  is 

exactly  what  we  expect  for  a  dielectric  medium  with  a  real 

p   p 
dielectric  constant  (l  -  cd  /co  ). 

We  conclude  this  section  by  obtaining  the  solution  in 

the  case  of  a  specular  reflection  boundary  condition.   As  is 

physically  reasonable,  but  by  no  means  necessary  mathematically, 

we  now  assume  that  f  (v)  is  even  and  we  look  for  a  solution 

o^  ' 

for  which  h  (v,t)  =  h~(-v,t)  (we  no  longer  restrict  h  (v,t) 
to  the  time  dependence  e~    ,  thus  H+(z  )/(p+icD)  now  becomes 
H(z,p).   We  start  from  (25)  perform  two  changes  in  the  paths 
of  integration.   First,  we  treat  the  term  proportional  to  E 


o 


exactly  as  above.   Second,  we  treat  the  two  integrals 
containing  the  terms  H(t,p)e~'  '^^'P\    H(t,p)e"^  (t,p) 
separately  and  use  the  contours  of  integration  of  Fig.  3, 
and  after  rewriting  the  result  we  find 
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pE    (en(z,p)-r(0,p)    _    ^) 

(32)      L(z,p)   =   H(z,p)   +  K(z,p)   =  — ^ 

z(p+lco) 


+ 


nz,p) 

27ri 


/^ 


(di_)g-r(t,p)(H+(^^p)_H-(t,p)) 


\7 


+ 


n(z,p) 

27ri 


r 


<j 


dt        ^G-(t,    p/t)e-n"(t. 


Xt^^ 


P, 


o 


-   G+(t,    p/t)e-r'''(t,p)  I 


The  condition  of  specular  reflection  Is  just  L(z,p)  =  L(-z,p), 
and  thus  we  must  determine   H  (t,p)  -  H  (t,p)  =  -27rith  (t,p) 
so  that  the  right  hand  side  of  (32)  Is  an  even  function  of  z. 
We  may  reduce  this  problem  to  a  very  trivial  Wlener-Hopf 
problem.   We  take  the  condition  L(z)  =  L(-z)  and  multiply  by 
e-r(z,p)-r(-z.P)  3o  that 

(|)  f^o\Cri-z,T^)-rio,v)  _  3-r(z,p)  -  n(-z.p)^ 


+ 


(33) 


27rl 


o 
r 


(^)e-r'(t.P)(H+(t,p)-H-(t,p)) 


—  00 


+ 


00 

^{j^)fG~it,    p/t)e-P"(^^P)  -  G+(t,  p/t)e-^7  = 


t-z 


27  - 


-  (?) 


z '   \  p+ico 


n(z,p)-r(o,p)     -r(z,p)-r(-z,p) 


+ 


rr(z,p) 

27rl 


r 


(||^)e-r'(t,p)(H+(t,p)-H-(t,p)) 


+    Mffe^r^'^^^    p/t)e-r"(^'P) 


■,+ 


G    (t,    p/t)e 


-r 


+^ 


We  note  that  we  may  easily  split  (35)  Into  the  sum  of  two  sets 
of  terms;  the  first  consists  of  functions  analytic  except  for 
a  cut  on  the  positive  real  axis  and  the  second  contains 
functions  analytic  except  for  a  cut  on  the  negative  real  axis. 
Both  sets  of  functions  are  bounded  at  the  origin  and  go  to  zero 
at  Infinity;  hence  both  sets  are  Identically  zero.   Me  carry 
out  the  details  for  the  case  f  (v)  and  hence  P(z)  even  and 
G(z,s)  =  0.   Since  the  Initial  data  do  not  affect  the  solution 
for  large  time  we  do  not  make  a  serious  error  In  their 
omission.   We  use  the  Identity 


P(z,p)   +    P(-z,p) 


1 


2Trl 


r 

/dt    >,-, 


/p2+a32tV(t) 


^ 


J 

—00 


V 


P^+a3pt^F~(t 


or 


(54)      r(z,p)  +  r(-z,p)  =  iog(p2+zVp(z)), 


with  these  assumptions  and  {j)h) ,    (53)  becomes 
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E 


o    (£)(e-n(-z,p)_g-r(o,p))g-r(o,p) 


(p+ico)    ^z 


+ 


r(-z,p) 

2irl 


r 


\j 


(_M)e-r(t,p)(H+(t,p)_H-(t,p)) 


pE 
-^   o 

7rl(p+la)) 


o 

r 


(35) 


dt 


v^ 


^^^^/   p2+a.2tV(t) 


p2+a32tV(t) 


L^  r.-P(2^p)  «-r(o,p),  -P(0,p) 


z  ( p+lco ) 


( 


-e 


)e- 


+ 


,-r(z,p) 

27rl 


(_g)e-n(t,p)(H+(t,p)-H-(t,p)) 


vy 


+ 


pE 
_    o 

7rl(p+lcD) 


00 

r 


dt 
t(t-z) 


vj 


p2+a.2t2p+(t) 


P^+OD^tV(t) 


The  left  hand  (right  hand)  side  of  (35)  has  a  cut  on  the 
negative  (positive)  real  axis.  Is  bounded  at  the  origin  and 
tends  to  zero  at  Inflnityo   Hence  both  sides  of  (35)  are 
entire  and  vanish  at  Infinity  so  that  they  are  Identically 
zero.   After  straightforward  manipulations  we  find 

o 


(36) 


E 


^(^'P)  =  s(p+lco) 


1.4 

TTl 


dt 


1 


\J 


t(^t-p)  1p2^^2,2p+(^) 


P 


1 


p^+co^t^F  (t 
^    P 


1]  . 
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We  may  pass  to  the  limit  of  large  time  exactly  as  before  and 
we  then  easily  find  the  standard  result 


(37)   lim  (E(x,t)e+^"^t)  =^  -  ^°    ^° 


o 
r 

dt'e-^^A 


,   s  2    TTi 

t  ->  00  CD 

CO         -00 


t'     1  l-oOpt'^P  (t'co) 


1-CO^t  ^P^(t  O)) 


IV.   Asymptotic  Expansions  of  the  Steady  State  Solution 

In  this  section  we  consider  the  asymptotic  form  of  the 
electric  field  after  a  long  time 

(38)        ^U,t,o,)   =   e"^^^  lim  (e(x,  t  )e^^^  ) , 

t  -^  00 

where  the  Laplace  transform  in  space  of  the  limit  is  given  by 
(30).   We  have  obtained  this  steady  state  solution  by  passing 
to  the  limit  of  large  time  in  a  time  varying  solution.   However, 
we  may  now  easily  characterize  this  solution  uniquely  by  the 
condition  that  (£^(x,  t  ,co)  possess  an  analytic  continuation  into 
the  half -plane  ,  Im  cd  >  0.   Such  a  condition  leads  to  the  same 
Hllbert  problem  (23)  (with  G  ~  0,  and  the  factor  (p+ico)  omitted) 
and  ultimately  to  the  same  solution  (30).   Such  an  analyticity 
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property  Is  typical  of  causality  conditions  In  scattering 
theory. 

We  have  already  obtained  one  approximation  f  or  £  (x,  tjco) 
for  large  x. 


(31    )  lim  6(x,t,a.)   =  f- 2-TT   (^o    -  ^^' 

X    ->  00  (1-    60/03     ) 


where 


J  =   47rq 


00  u 

vh    (v)dv  +        vh~( V, t=oo)dA 


o 


but  we  would  like  sharper  estimates  analogous  to  the  results 
In  Section  II  for  the  Initial  value  problem.   Earlier  we  were 
able  to  derive  expansions  on  the  basis  of  two  distinct 
assumptions.   First,  we  were  able  to  adjust  the  contour  of 
Integration  In  (3)  so  that  we  could  take  uniformly  valid 
asymptotic  expansions  of  all  the  functions  In  question.   In 
order  that  the  error  terms  be  small  we  required  that  t  be  small. 
In  the  present  case  the  contour  of  Integration  for  the  Inversion 
of  the  Laplace  transform  of  (°(x,t,cD)  passes  through  s  =  00,  and 
so  the  only  expansions  that  we  can  hope  to  be  uniformly  valid 
are  those  valid  near  s  =  00.   We  could  Indeed  find  such 
expansions,  which  apply  for  small  values  of  x,  by  taking  a 
contour  of  Integration  far  from  the  origin  In  the  s-plane  and 
expanding  the  functions.   However  these  expansions  are  probably 
none  too  significant,  as  the  results  depend  very  sensitively  on 
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the  precise  behavior  of  the  distribution  functions  for  small 
velocity.   Further,  they  do  not  exhibit  any  plasma  oscillation 
phenomena  and  hence  we  omit  them.   Our  second  method  of 
obtaining  expansions  in  section  II  relied  on  the  observation 
that  the  integral  (3)  was  singular  at  kA  =  0,  and  hence  there 
must  be  approximations  possible  for  kA  small.   Here  the 
analogous  observation  is  that  (26),  (30)  and  (31')  are  singular 
at  CO  =  CO  and  hence  we  expect  interesting  results  for  1 


2/  2 


'P  ^  "  P 

small.   We  treat  this  limiting  case  in  this  section  in  two  parts, 

for  we  must  examine  separately  1  -  <^  /co   small  and  positive  and 

then  small  and  negative.   We  obtain  results  relevant  for  a 

large  enough  range  of  values  of  x  to  observe  Landau  damping 

of  plasma  oscillations.   As  a  matter  of  convenience  we  restrict 

ourselves  to  type  I  distribution  functions. 

Our  first  task  is  to  find  appropriate  representations  for 

P(z,-ico)  for  all  real  co.   We  start  from  (26a),  which  provides 

P    P 

the  analytic  continuation  of  [""(z^-ico)  for  all  co,  co  7^  co  ,  and 

we  wish  to  undo  the  transformation  that  led  to  (26a)  and  return 
to  a  form  analogous  to  (24).   With  p  =  -ico,  (26a)  becomes 


(26'a)  r(z,-icD) 


1 


27ri 


dt 


(t  -  z/cu) 


log 


(l-co'^t  "^P(cot  ))co^ 

2   2 
CO  -co 

p 


+  l0£ 


O  O  — 

'  CO"  -co  z  P '  (  z  ) 

_^ 

2   2 
a.  -cOp 
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where  we  have  taken  Im  z  <  0  (and  we  assume  throughout  that 
03  >  o).   The  reader  may  readily  verify  that  whenever  P(z,-lco) 
Is  used  In  (30)  It  corresponds  to  values  of  z  in  the  lower 
half-plane  or  on  Its  boundary,  and  whenever  there  Is  an 
ambiguity  for  real  o),  the  result  is  the  limiting  value  for 
CO  in  the  upper  half-plane. 

We  wish  to  take  the  section  of  the  integral  from  minus 
Infinity  to  zero  and  deform  it  through  the  lower  half-plane 
into  the  path  from  zero  to  plus  infinity.   Besides  the  pole 
in  the  lower  half-plane  at  t  =  z/co  there  may  be  a  branch  cut 
at  a  zero  of  the  logarithm.   For  a  stable  plasma  it  is  readily 
shown  that  in  the  lower  (or  upper)  half -plane  t  F(t)  assumes 
all  real  values  from  zero  to  one  exactly  once.--^  We  denote 

by  t~  the  zero  of  (co  -co  t  P(t))  in  the  upper  or  lower  half-plane, 

2    2 
CD  <  CO  .   To  accomplish  the  transformation  described  above  we 

use  the  contour  of  integration  of  Pig.  4  and  we  find 


5.   Prom  the  asymptotic  expansion  it  is  clear  that  there  is 
a  curve  near  the  imaginary  axis  at  infinity  on  which 
t  P-(t)  is  real  and  positive.   It  is  easy  to  show  that 
the  curves  enter  the  origin  and  there  are  no  other  such 
curves.   Hence  t  P~(t)  is  monotone  on  the  curve  and 
takes  all  real  values  from  O(t  =  0)tol(t  =  oo). 
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(39) 


P(z,-ia))  = 


(ffi)  log 


'cD2-a)2tV(t) 


CD^-a)^t^P  (t) 


+ 


log(- 


z-t 


CDn 


0  <  CD  <  CO 


0 


0^    >    0)  J 


and  In  (39)  the  logarithm  tends  to  zero  as  t  tends  to  Infinity, 
Hence   P(z,-lto)  Is  analytic  In  z  except  for  a  cut  on  the 
positive  real  axis  and  possibly  a  branch  point  In  the  lower 
half -plane.   We  see  below  that  P(z,-la))  Is  In  fact  bounded 
for  small  z  even  though  there  Is  an  apparent  singularity  at 
z  =  0.   With  this  Information  we  may  modify  the  contour  of 
Integration  In  (30 )  and  finally  write 


(^0) 


where 


e        dsx{sj, 


Re   s  >   0 


(41)      I(s) 


r^ 


ICD  .       -i 

-     ICO  J 


s(l+%F(-f^) 


-n(o,-io3) 


o 

r 


+ 


27rl 


dt'(H''(a3t'}-H-(:a3t'))e-^(^^'-^^)  1 


—  00 


(st    +1) 
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We  readily  Infer  that  l(s)  is  analytic  in  s  except  for  a  cut 
on  the  positive  imaginary  axis,  a  pole  at  s  =  -ico/t"'" 


CO 


(0  <  o)  <  od   )  ,    and  a  pole  at  s  =  0.   We  deform  the  contour 
according  to  Pig.  5  to  obtain 


(E   -  ^)e-^^t 
(^2)  £(x,t)    =        °      ^%      ^ 

(1     -    03^0.2^ 


100 


-icot 


+ 


27ri 


e^''ds(l+(s)    -   I"(s)) 


o 


+ 


lim 


-ioo/t 


CD 


i^  +  ^   )l(s)      e 


OicD(t+  x/t    ) 


CD' 


t 


CO 


0 


0  <   CO  <  CO. 


P 


,  CO    >    CO 


+  , 


We  note  that  if  both  f  (v)  and  h  (v)  are  type  I  distribution 


o 


.+ 


functions  then  (l  (s)  -  I~(s))  =  Qi{s    )    for  s  small  for  any  n 
and  l(s)  =  O^i")  ^°^  ^  large.   Hence  for  cd  fixed  and  different 
from  CO  the  expansion  of  E(x,t)  for  large  x  is  just  (31') 
with  an  error  term  0(1/^  )  obtained  by  repeated  integration 
by  parts.   Again,  we  believe  this  to  be  the  only  expansion 
possible  for  large  x  which  does  not  depend  on  the  detailed 

properties  of  the  distribution  functions  and  boundary  condl- 

2   2 
tlons.   In  order  to  obtain  expansions  for  1  -  <^^<^  small  we 

must  return  to  (39)  and  obtain  asymptotic  expansions  for  large 

2   2 
z  uniform  in  (1  -  to^<^  )• 
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p    p 

We  consider  first  the  simpler  case  of  (l  -  cd  /co  )  small 

+  ,  . 
negative.   Prom  the  asymptotic  expansions  of  P~(t)  we  see 

that  CD  -CD  t  P"(t)  Is  approximately  real  and  negative  for  t 

large  and  real  so  that  the  logarithm  In  the  Integrand  In  (59) 

Is  uniformly  0(l/t  )  for  any  n.   Thus  the  asymptotic  expansion 

gives  no  difficulties  and  we  find  easily 

(^3)     r(z,-la3)  =  21  ^  +0(1A^^-)  +  log(^) 

n=l   z 


where  the  coefficients  and  error  estimate  are  uniformly  bounded 

p   p 
In  (l  -  CD  /oD  ).   The  standard  methods  also  show  that  the 

Integral  In  (59)  is  log  z  +  0(  1 )  fo^  small  z,    so  that 

P(zj-Icd)  1s(3(1)  for  small  z.   In  order  to  complete  our 

discussion  of  the  expansion  of  the  electric  field  in  this  case 

+ 
we  need  an  approximation  to  t-  and  we  find  easily  that 


For  this  case  we  almost  have  the  complete  answer  in  ( 42) 
since  we  have  an  explicit  answer  plus  an  Integral  which  we 
wish  to  show  small.   We  divide  the  integral  into  two 
Integrals,  one  from  0  to  IRcd  and  then  other  from  IFkjo  to  loo. 
For  the  first  Integral  we  take  the  asymptotic  expansion  to 

th 

N   order  in  s  (which  is  identically  zero)  plus  an  error 


56 


term.   In  the  second  Integral  we  Integrate  by  parts  repeatedly 

n  times  and  evaluate  the  end  point  terms  approximately  by  the 

t  h 
asymptotic  expansion  to  N   order.   If  we  use  the  expansion  (44) 

explicitly  only  to  lowest  order  then  we  find 


(^5)    8U,t)   =  ^^°  "  ^%^      (l-exp, 


tax 


(1-  mW)^ 


+ 


""  ^r  X-"  (1-  o?/is?)~'^ 


,T3m+l//T    2/  2>i,o2\n-m+l 
m=o    |R    ((1-  o)  /CD  j+R  ) 


In  Eq.  (45)  we  have  by  no  means  the  sharpest  error  estimate 
possible,  but  it  suffices  to  show.  Just  as  in  Section  II,  that 
the  exponentially  damped  component  of  the  electric  field  is  in 
fact  observable  until  it  damps  significantly  and  then  the  non- 
constant  part  of  the  electric  field  disperses  in  an  irregular 
manner.   We  omit  the  details  on  the  choices  of  R,  N,  and  n 
since  they  are  essentially  the  same  as  the  discussion  following 
in  Eq.  (l4).   The  exponential  damping  observable  in  this  case 
is  not  "Landau  damping"  since  it  is  not  associated  with  plasma 
oscillations  and  since  it  would  be  present  in  a  macroscopic 
theory  as  well  as  in  the  present  kinetic  approach. 
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We  now  consider  the  case  1  -  <^„/co   small  and  positive. 
The  expansion  of  P(z,-lcu)  for  large  z  is  somewhat  more 
complicated.   Since  (F"''(t))   =  P~(t)  we  Infer  that 


(39')      n(z,-ico)  =  i 


00 

r 

dt  „  „/  2   2^2t^+, 


,         arg(co^-co^t^F'^(t)), 
t-z   ^^  p     \  /  / 


o 

Prom  the  asymptotic  expansion  of  F  (t)  we  Infer  that 
Re(cD  -co  t  F  (t))  changes  from  negative  to  positive  as  t 

p     p   P   _|_ 

Increases  and  is  large.   In  fact  Re(a)  -(^„t  F  (t))  changes 

sign  at  least  one  other  time  since  the  value  is  positive  at  t 

equals  zero,  but  we  are  concerned  with  the  changes  in  sign 

+ 
in  the  asymptotic  region  only.   We  denote  hy  t   the  positive 

Q      P  — 

(negative)  value  of  t  approximately  equal  to  ±(a)  /co   -  l)^ 
for  which  Re(co^-co^t^P^(  t ) )  =  0.   Since  f^lv)  is  a  type  I 
distribution  function  it  follows  that  in  any  but  an  exceedingly 
narrow  interval  surrounding  t 

^       CO 

arg(m''-mpt2F+(t))  = 

-  +  0  {71)  t  <  *: 

The  width  of  the  interval  in  which  the  above  estimate  fails  is 
less  than  (^(l/t  )  for  any  N.   Thus  if  we  add  and  subtract 
TrT](t  -t)  from  (39')  we  readily  obtain  the  expansion 
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n=l  z       z   '      ^  ' 


valid  for  z  bounded  away  from  the  positive  real  axis  and 
again  the  coefficients  are  uniformly  bounded  in  (l  -  cd  /cd  ). 
We  note  the  formal  similarity  to  (43)  for  the  case  cd  <  cd  . 
Again  the  standard  arguments  show  that  P(z,-ia))  is  Q){l)   for 
small  z. 

We  are  now  prepared  to  obtain  the  asymptotic  expansion 
of  the  electric  field.   For  the  values  of  x  which  are  not  too 
large  we  use  the  contour  of  Pig.  6,  where  the  circle  is 
centered  at  s  =  -ioo/t"  and  is  of  radius  p,    and 

'  OD  ^ 


g(x.t).(^-l^)C 


-ICDt 


{l-a>2/(«2)    2Tri 


e^^  dsl(s). 


We  now  follow  the  same  procedure  as  before.   We  split  the 
integral  into  two  parts,  an  integral  from  0  to  IRcd,  R  >  p, 
and  an  integral  from  IRco  to  infinity.   On  the  first  part  of 
the  contour  we  use  the  asymptotic  expansions  to  N   order 
terms.   On  the  second  part  of  the  contour  we  integrate  by 
parts  n  times  and  use  the  asymptotic  expansions  to  estimate 
the  error  terms.   If  we  carry  the  explicit  answer  only  to 
the  lowest  order  non-trivial  terms  we  find 
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(47) 


(f(x,t) 


(E   -    J/lco) 


(1-  coVco^) 
\  p/       / 


2-^(  l^exp/loox   /^^^-^    lU"^"'^ 


m=o 


From  (47)  we  easily  Infer  that  we  may  follow  the  undamped  plasma 
oscillation  for  a  distance  of  the  order  of  a  reciprocal  power 
of  (cD  /co   -  l)»   Again  we  need  a  more  sensitive  expansion  to 
obtain  any  damping.   We  use  the  contour  of  Pig.  6  again  except 
that  we  move  the  semi-circle  on  the  right  half-plane  onto  the 
imaginary  axis.   We  use  a  Taylor  series  plus  remainder  on  the 
sections  of  the  integral  from  -Ico/t   -  ip  to  -ico/t   +  ip,  and 
integrate  by  parts  n  times  on  the  remainder  of  the  contour. 
We  move  the  section  of  the  contour  on  which  the  Taylor  series 
was  applied  onto  the  semi-circle  in  the  left  half-plane, 
getting  a  contribution  from  the  pole  near  s  =  -ico/t",  and  we 
ln1:egrate  by  parts  n.  times.   We  thus  obtain 
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(^8)   £(x,t) 


(E^-  J/lco) 
(1-  co2/a)2) 


1-exp  xo)  (1 


(^)  ('^   -1 


-icot 


n-m+1   n-m+l| 


where 


(49) 


5  =  -  ^ 


IT  /       O^  ' 

"5u 


U 


6a>5 


r,  2  r  2/  2    r 

2CD„   /O)  /CD.,  -1 


P 


P 


If  In  (48)  we  take  -^  =  ■ifcoV'^  -l)^  ^  then  the  result  Is 
essentially  the  same  (l6)  with  kA  replaced  by  (cu  /co  -l)^  , 
and  again  it  is  easy  to  show  that  it  is  possible  to  observe 
the  damping  of  the  plasma  oscillations  long  enough  to  see 
their  exponential  damping,  but  after  the  net  damping  is 
large  there  is  only  a  slow  dispersion  of  the  electric  field 
to  its  equilibrium  value  (31  )» 
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Figurel  A  Contour  of  Integration 


Figure  2  A  Contour  of  Integration 
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Figure  4  A  Contour  of  Integration 


Figure  5  A  Contour   of   Integrotion 


sz 


Figure  6  A   Contour  of  Integration 
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